Quasi-one-dimensional Bose gases with large scattering length 
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Bose gases confined in iiigiily-elongated iiarmonic traps are investigated over a wide range of 
interaction strengths using quantum Monte Carlo tecliniques. We find tfiat tfie properties of a Bose 
gas under tiglit transverse confinement are well reproduced by a Id model Hamiltonian with contact 
interactions. We point out the existence of a unitary regime, where the properties of the quasi- Id 
Bose gas become independent of the actual value of the 3d scattering length a^d- In this unitary 
regime, the energy of the system is well described by a hard rod equation of state. We investigate 
the stability of quasi-ld Bose gases with positive and negative a^d- 
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In recent years the study of quasi-ld quantum Bose 
gases has attracted a great deal of interest. Intriguing 
properties of quasi-ld gases, such as the exact mapping 
between interacting bosons and non-interacting fermions, 
have been predicted P, ■ ^ bosonic gas that behaves 
as if it consisted of spinless fermions, a so-called Tonks- 
Girardeau (TG) gas, cannot be described within mean- 
field theory since it exhibits strong correlations; instead, 
a many-body framework is called for. While experimen- 
tal evidence of quasi-ld behavior has been reported for 
bosonic atomic gases under highly-elongated harmonic 
confinement 0], TG gases have not been observed yet. 
It has been suggested, however, that TG gases can be 
realized experimentally for either low atomic densities or 
strong atom-atom interaction strengths. The 3d s-wave 
scattering length asd, and hence the strength of atom- 
atom interactions, can be tuned to essentially any value, 
including zero and ±oo, by utilizing a magnetic atom- 
atom Feshbach resonance 0, . 

Utilizing a two-body Feshbach resonance, 3d degener- 
ate gases with large scattering length 03^ have been stud- 
ied experimentally and theoretically. For 03^ ±00, it is 
predicted that the behavior of the strongly-correlated gas 
is independent of 03^ 0. For homogeneous 3d Bose gases, 
this unitary regime can most likely not be reached ex- 
perimentally since three-body recombination is expected 
to set in when 03^ becomes comparable to the average 
interparticle distance. Three-body recombination leads 
to cluster formation, and hence "destroys" the gas-like 
state. The situation is different for Fermi gases, for which 
the unitary regime has already been achieved experimen- 
tally 1^. In this case, the Fermi pressure stabilizes the 
system even for large \a3d\- It has been predicted that 
three-body recombination processes are suppressed for 
strongly interacting Id Bose gases These studies 

raise the question whether a highly-elongated inhomoge- 
neous Bose gas, that is, an inhomogencous quasi-ld Bose 
gas, is stable as a^d — > ±00. 



This Letter investigates the properties of a quasi-ld 
Bose gas at zero temperature over a wide range of in- 
teraction strengths within a microscopic, highly accurate 
many-body framework. We find that the system i) is well 
described by a Id model Hamiltonian with contact inter- 
actions and renormalized coupling constant for any 
value of the 3d scattering length a^d] ii) behaves like a 
TG gas for a critical positive value of 03^; iii) reaches a 
unitary regime for large values of |a3d|, where the proper- 
ties of the quasi-ld Bose gas become independent of the 
actual value of 03^ and are similar to those of a hard-rod 
gas; and iv) becomes unstable against cluster formation 
for a critical negative value of a^d- 

Our study is based on the 3d Hamiltonian H^di 
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which describes N spin-polarized mass m bosons un- 
der highly-elongated confinement with ujz = Awp, where 
A <C 1. The coordinates pi — \J ^1 ^ yf and Zi denote, 
respectively, the radial and axial position of the «-th par- 



ticle, rij — 



denotes the interparticle distance 



between atom i and j, and V{r) denotes the two-body 
interatomic potential. We consider two different poten- 
tials: i) a purely repulsive hard-sphere (HS) potential, 
V-^^{r) = 00 for r < 03^, and zero otherwise; and ii) a 
more realistic short-range (SR) potential, which can sup- 
port two-body bound states, V^^'{r) = — Vb/cosh^(r/ro), 
with well depth Vq. For V^^, the s-wave scattering 
length 03^ coincides with the range of the potential. For 
V^^, in contrast, rg determines the range of the poten- 
tial, while the scattering length 03^ is a function of rg 
and Vq. In our calculations, ro is fixed at a value much 
smaller than the transverse oscillator length, vq — O.lOp, 
where Op — y^h/mujp. To simulate the behavior of a^d 
near a field-dependent Feshbach resonance, we vary the 
well depth Vq and consequently the scattering length 03^. 
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FIG. 1: gid [dashed line, Eq. @] and aid [solid line, Eq. 
as a function of astj. A vertical arrow indicates the value of 
asd where gid diverges, a^d/cLp = 0.9684. Horizontal arrows 
indicate the asymptotic values of gid and aid, respectively, as 
asd ±00 {gid = —1.9368aphu!p and aid 
asd as a function of the well depth Vo for V^^{r). 



Importantly, 03^ diverges for particular values of Vq (the 
inset of Fig. Q shows one such divergence). At each di- 
vergence, a new two-body s-wave bound state is pulled 
in. Our numerical calculations (see below) are performed 
for the well depths Vq shown in the inset of Fig. ^ 

We consider situations where the bosonic gas described 
by Eq. is in the Id regime for any value of the 3d 
scattering length 03^, which implies A^A <C 1. If the 
range of V{r) is much smaller than Cp, it is predicted (21 
that the properties of the 3d gas are well described by 
the Id contact-interaction Hamiltonian Hid, 
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where gid is an effective coupling constant. The limit 
gid — > corresponds to the weakly-interacting mean- 
field regime, while gid — > 00 corresponds to the strongly- 
interacting TG regime. For positive gid and vanishing 
confinement {ojp = ujz = 0), Eq. reduces to the 
Lieb-Liniger (LL) Hamiltonian 9], whose gas- like prop- 
erties have been studied in detail. For negative gid and 
uJn = ujz = 0, Eq. (01 supports cluster-like bound states 
[lil; little is known about gas-hke states in this case. 

Below, we solve the many-body Schrodinger equation 
for the Id Hamiltonian with confinement, Eq. for pos- 
itive and negative gid, and relate our results to those for 
the 3d Hamiltonian, Eq. Considering two bosons in 
a highly-elongated geometry that interact through a reg- 
ularized zero-range pseudo-potential, Olshanii [2| shows 
that the effective Id coupling constant gid can be ex- 
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where A = |C(1/2)|/V^ = 1.0326. Alternatively, gid can 
be expressed through the effective Id scattering length 
aid, gid = -2Ti' / {maid) 2], where 



aid 



A 



asd 



(4) 



For 03^ > (but a^d ^ o-p), gid approaches the un- 
renormalized coupling constant, 17°^ = '^'h^o-^d/ {'ma^p), 
which is obtained by averaging the 3d coupling constant 
g'id — 'inh'^ a^d / m over the transverse oscillator ground 
state (see, e.g., Ref. 0). For these values of asd, the Id 
Hamiltonian given in Eq. |(5J) with gid replaced by g^d, 
describes a quasi-ld system accurately 



For asd 



> 



repl; 



and for negative a^d, in contrast, the con- 



finement induced renormalization becomes important, 
and the effective Id coupling constant gid and scatter- 
ing length aid, Eqs. Q and Q), have to be used. Fig- 
ure n shows gid (dashed line) and aid (solid line) as a 
function of a^d- At the critical value = 0.9684ap (in- 
dicated by a vertical arrow in Fig.^, 171^ diverges while 
aid goes through zero. At the 3d resonance, that is for 
asd ±00, gid and aid each reach an asymptotic value 
{gid — —l.QSGSapfiujp and aid 



1.0326ap, respectively. 



indicated by horizontal arrows in Fig. Tuning 03^ 
to large values hence allows a unitary quasi-ld regime, 
where 171^ and aid are independent of 03^, to be entered. 

We solve the Schrodinger equation for the many-body 
3d Hamiltonian, Eq. (Q, numerically using the diffusion 
quantum Monte Carlo (DMC) technique. Our interest 
is in the lowest gas- like many-body state. For V^^ , the 
lowest gas-like state of H^d coincides with the many-body 
ground state; to describe this state we can hence use the 
"standard" DMC technique (see, e.g., For V'^^, 

in contrast, the many-body ground state is a cluster-like 
bound state. To describe the lowest-lying gas-hke state, 
i.e., an excited state, we hence use the fixed- node DMC 
(FN-DMC) method 131. For a given many-body nodal 
surface, the FN-DMC method allows the Schrodinger 
equation to be solved for approximate excited states. For 
N = 2, we obtain the exact nodal surface by direct di- 
agonalization. Assuming that the scattering properties 
between each atom pair are unaltered by the presence 
of other atoms, we then use the N — 2 nodal surface 
to construct the many-body nodal surface. For dilute 
quasi-ld gases we expect that the FN-DMC approach as 
implemented here results in highly accurate many-body 
energies. 

Figure 2 shows the resulting 3d energy per particle, 
E/N — hiOp, as a function of a^d for = 5 under quasi- 
ld confinement, A = 0.01, for the hard-sphere and the 
short-range two-body potential, V^^ (diamonds) and 
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FIG. 2: 3d DMC energy per particle calculated using V^^ 
(diamonds) and V^^ (asterisks), respectively, together with 
Id DMC energy per particle calculated using g^d [squares, 
Eq. Q] and g^^ (pluses), respectively, as a function of [(a) 
linear scale; (b) logarithmic scale] for N — 5 and A — 0.01. 
The statistical uncertainty of the DMC energies is smaller 
than the symbol size. Dotted and solid lines show the Id 
energy per particle calculated within the LDA for gl^ (using 
the LL equation of state) and for gid, Eq. Q (using the LL 
equation of state for gid > 0, and the HR equation of state 
for gxd < 0), respectively. A dashed horizontal line indicates 
the TG energy, and a vertical arrow the position where gid, 
Eq. Q, diverges. 



V^^ (asterisks), respectively. For small azd/o-p, the en- 
ergies for these two two-body potentials agree within the 
statistical uncertainty. For 03^ ^ Op, however, clear dis- 
crepancies are visible. The DMC energies for V^^ cross 
the TG energy per particle (indicated by a dashed hor- 
izontal line), E /N — TiLOp = huipXN/2, very close to the 
value Og^j = 0.9684ap [indicated by a vertical arrow in 
Fig- nib)], while the energies for V^^ cross the TG en- 
ergy per particle at a notably smaller value of 03^. 

To compare our results obtained for the 3d Hamilto- 
nian, -ffsd, with those for the Id Hamiltonian, Hid, we 
also solve the Schrodinger equation for Hid, Eq. For 
positive values of the coupling constant 171^ we calculate 
the many-body ground state energy by the exact DMC 
method. For gid < 0, however, the Id Hamiltonian sup- 
ports many-body bound states, and, as in the 3d case, 
we use the FN-DMC method to describe the lowest-lying 
gas-like state. For N = 2 and gid < 0, the first excited 
state of the Schrodinger equation for Hid, Eq. has a 
node at 2:12 = aid- To solve the many-body Schrodinger 
equation for Hid with negative gid for the lowest gas-like 



state by the FN-DMC technique, we parameterize our 
many-body nodal surface in terms of aid- This many- 
body nodal surface is expected to be good if the density 
of the gas is low. 

In addition to the 3d energy per particle. Fig. 2 shows 
the resulting Id energy per particle obtained by solving 
the Schrodinger equation for Hid, Eq. Q, for the renor- 
malized coupling constant gid [squares, Eq. ||2J)], and the 
unrenormalized coupling constant g°^ (pluses), respec- 
tively. The Id energies calculated using the two different 
coupling constants agree well for small 03^, while clear 
discrepancies become apparent for larger 03^. Impor- 
tantly, the Id energies calculated using the renormalized 
Id coupling constant gid agree well with the 3d energies 
calculated using the short-range potential V^^ (aster- 
isks) up to very large values of the 3d scattering length 
03^, and also for negative 03^. In contrast, at large a^d 
the Id energies deviate clearly from the 3d energies calcu- 
lated using the hard-sphere potential V^^ (diamonds). 
We conclude that the renormalization of the effective Id 
coupling constant gid and the Id scattering length aid 
are crucial to reproduce the results of the 3d Hamilto- 
nian Hjd when 03^ Op and when a^d negative. Small 
deviations between the Id energies calculated using the 
renormalized Id coupling constant gid and the 3d ener- 
gies calculated using the short-range potential V^-'^ re- 
main; we attribute these to the finite range of V^^. 

If the size of the cloud is much larger than the har- 
monic oscillator length a^, where az = yjfi/mujz, it has 
been shown that the properties of the Id LL Hamiltonian 
Hid, gid > 0, are well described by a simple equation of 
state using the local density approximation (LDA) [l^. 
For gid < 0, we instead apply the equation of state for Id 
hard-rods (HRs). Recall that the many-body nodal sur- 
face of the lowest-lying gas-like state of Hid with gid < 
is well parameterized by aid- For > aid, the corre- 
sponding wave function coincides with that of N Id hard 
rods of size aid- For small values of the Id gas parame- 
ter, nidaid ^ 1, where nid denotes the linear density, we 
hence expect that the lowest-lying gas-like state of the Id 
many-body Hamiltonian with gid < is well described 
by a system of hard rods of size aid- The exact energy 
per particle of the uniform hard rod system is given by 
E/N = {Tr'^h^nlJ6m)/il - nidaidf 0. For trapped 
systems with iVA ^ 1 we obtain the expansion 



iVA 



E/N — hujp = hujp — 



128V2 
457r2 ' 



-aid 



.(5) 



The first term corresponds to the energy per particle in 
the TG regime. Similarly, the linear density in the center 
of the cloud is to lowest order given by the TG result, 
Jiid = V2NX/ {trap)- In the unitary limit, that is for 
aid/ap = 1.0326, expression © becomes independent of 
a^d, and depends only on NX. 

Lines in Fig. [21 show the resulting Id energies per par- 
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tide for the LL equation of state {gid > 0) as well as 
the HR equation of state {gid < 0) calculated within 
the LDA. Remarkably, the LDA energies nearly coin- 
cide with the Id many-body DMC energies (pluses and 
squares, respectively); finite-size effects play a role only 
for asd <C a p. Our calculations establish for the first 
time that a simple treatment, i.e., a HR equation of 
state treated within the LDA, describes trapped quasi-ld 
gases with negative coupling constant gid well over a wide 
range of the 3d scattering length a^d- For a^d 0, that 
is, for large aid, the HR equation of state using the LDA, 
cannot properly describe trapped quasi-ld gases, which 
are expected to become unstable against formation of 
cluster-like many-body bound states for aid ~ ^/nid- We 
hence investigate the regime with negative a^d in more 
detail within a many-body framework. 

By comparing with results for the 3d Hamiltonian, 
Eq. we have shown above that the Id Hamiltonian 
Hid, Eq. provides an excellent description of quasi- 
ld gases. We hence base our stability analysis of quasi- 
ld gases with large effective Id scattering length aid on 
Hid- We solve the many-body Schrodinger equation for 
Hid using the variational quantum Monte Carlo (VMC) 
method. Our variational A^-particle Bijl-Jastrow-type 
wave function consists of one- and two-body terms. The 
one-body terms are written as a function of a single vari- 
ational parameter a, which determines the size of the 
atomic gas. The two-body term is parameterized by aid, 
and explicitly accounts for correlations. 

Figure |21 shows the resulting VMC energy per parti- 
cle for N = 5 and A = 0.01 as a function of the vari- 
ational parameter a (Gaussian width) for four different 
aid- For aid/ap — 1 and 2, Fig.|21shows a local minimum 
at amin ~ ciz- The minimum VMC energy nearly coin- 
cides with the essentially exact DMC energy, which sug- 
gests that our variational wave function provides a highly 
accurate description of the quasi-ld many-body system. 
The energy barrier at a ~ 0.20^ decreases with increas- 
ing aid, and disappears for aid/ap « 3. We interpret this 
vanishing of the energy barrier as an indication of insta- 
bility jl5|]. For small aid, the energy barrier separates the 
gas-like state from cluster-like bound states. For larger 
aid, this energy barrier disappears and the gas- like state 
becomes unstable against cluster formation. 

We additionally performed variational calculations for 
larger TV and different A. We find that the onset of in- 
stability of quasi-ld Bose gases can be described by the 
product of the Id scattering length aid and the linear 
density at the trap center nid- To be specific, our many- 
body calculations suggest that a quasi-ld gas is stable for 
aidiT-id ^ 0.35, and becomes unstable for ai^n-id ?J 0.35. 
Our analysis suggests that the quasi-ld unitary regime 
can be reached experimentally. By tuning the 3d scat- 
tering length, it is further possible to investigate the on- 
set of instability. By reducing A one should be able to 
stabilize relatively large quasi-ld systems. 
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FIG. 3: VMC energy per particle as a function of the varia- 
tional parameter a for A*' = 5, A = 0.01, and aid/flp = 1.0326 
(pluses), 2 (asterisks), 3 (diamonds) and 4 (triangles). 



In conclusion, we investigated the energetics of a Bose 
gas under highly-elongated harmonic confinement over a 
wide range of the 3d scattering length. We find that the 
quasi-ld gas can be described by a many-body Id model 
Hamiltonian with zero-range interactions and renormal- 
ized coupling constant. For a^d ±cx), the quasi-ld gas 
enters a unitary regime, where all properties of the sys- 
tem are independent of a^d- In the vicinity of the unitary 
regime, the quasi-ld system behaves like a gas of HRs. 
For negative 03^, quasi-ld gases become unstable against 
cluster formation for a critical value of aidUid- 
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